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Abstract 

We present a calculation of the renormalized HQET Lagrangian at order 
C(l/m|) in the one particle sector. The anomalous dimensions of local 
operators and time ordered products of dimension 7 contributing at this 
order are calculated in the one loop approximation. We show that a 
careful treatment of the time ordered products is necessary to arrive at a 
gauge independent renormalized lagrangian. Our result sets the stage for 
an investigation of reparametrization invariance at 0(1/?tiq). 
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1 Introduction 



Heavy Quark Effective Field Theory (HQET) Q has been established as the 
theoretical tool of choice for the description of mesons and baryons containing 
one heavy quark M . This derives from the fact that it is a systematic expansion 
in inverse powers of the heavy quark mass with well defined and calculable 
coefficients. Furthermore, its realization of the spin and flavor symmetry of the 
low energy theory is a phcnomcnologically powerful tool. 

The 1 /uiQ expansion has already been applied sucessfully to phenomenologi- 
cal problems such as the determination of V c b- This involves a matrix element of 
the left handed current for the b — > c transition between two heavy hadrons. In 
order to match the experimental precision corrections (1/mg as well as radiative 
corrections) are indispensable. 

The present paper considers partial aspects of the 1 /toq corrections, which 
eventually are needed for phenomenological applications. Here the lagrangian 
of QCD is expanded up to order and the one- loop renormalization of 

the corresponding operators is calculated. Results for the lower order terms 
H ||, |[ ^, 0, |§] as well as a matching calculation for terms of order 1 /itiq 
are already known, but the full renormalization has not yet been studied. 

The lagrangian itself is also of theoretical interest, since on one side it is 
relevant for the 1/mg expansion of the heavy hadron mass, and on the other 
side allows a study of reparametrization invariance. 

To keep things as simple as possible we restrict ourselves to the one particle 
sector of HQET, disregarding four fermion operators as well as pure gluonic 
operators. The calculation of the anomalous dimensions of the corresponding 
operator basis sets the stage for an extensive study of reparametrization invari- 
ance at C(l/mg), which will be subject of a sequel to this paper. 

This note is organized as follows: in section ^|we introduce our operator basis 
and discuss its reduction to a set of linearly independent operators. Our result 
for the anomalous dimensions and some checks of its consistency are presented 
in section ||. In section ^ we give insight into some aspects of the calculation. Fi- 
nally, we present the logarithmic contributions to the short distance coefficients 
of the effective lagrangian in section |^ and conclusions in section ^. 



2 Operator basis 

The effective HQET lagrangian is defined by a systematic expansion of QCD in 
inverse powers of the heavy quark mass 

£hqet - h v (ivD)h v + E C i n) °i n) ■ (!) 

To lowest order only one operator, h v (ivD)h v , contributes, which is indepen- 
dent of the flavor and the spin of the heavy quark, resulting in the well known 
spin-flavor symmetry of HQET. However, these symmetries are valid only in the 
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limit rriQ — > oo and broken by local operators of ) of dimension 5 and higher. 
The coefficients C\ n ^ are short distance corrections, which compensate for the 
modified UV behavior of the effective theory. They are determined perturba- 
tively by matching HQET to QCD order by order in the strong coupling a s and 
the inverse heavy quark mass supplemented by a renormalization group analysis 
of the respective operators. 



We start with a review of the local operators appearing at 0(l/mg) and 
C(l/mg), respectively. At 0(l/rriQ) we choose the conventional basis 



0[ x) = h v (iDfh v 



O, 



(i) _ .9 



(2) 



Of ) = h v {ivDfh v 



and at 0(l/m Q ) 



(2) 

of) 



hyiD^ivD^D^hy 

Ily(ivD)(iD) 2 hy 



Of) = hyicr^iD^ivD^Dxhy 



o 



h v (iD) 2 (ivD)h v 



Oj' J = h v {ivDfhy 



(2) _ 
5 — 

Of) = h v i(j^ x iD^iD x {ivD)h v 



(2) 
4 

Of) = h v {ivD)i(j^ x iD^iD x h v 



(3) 



The definition of the covariant derivative iD = id + g s T a A a and field strength 
tensor F® v T a = —i/ g s [iD^,iD u ] follows usual conventions. Note, that we re- 
frain from skipping operators which vanish by the heavy quark equation of 
motion (EOM) (ivD)h v = 0. This will proove important, since once inserted 
into time ordered products at 0(l/rag), such operators contribute to physical 
operators, as will be explained in detail below. 

At 0(1/ tUq) we find 13 local operators contributing to physical matrixele- 
ments: 



Of 
Of 

of 
of 
of 



o 
o 



hyiD^ivDfiD^hy 

hyiD^iDfiD^hy 
hyiO^iD^ivDfiDvhy 

h v ia^ v iD p iD^iD v iD p h v 
hyia^iD^iDfiDvhy 
h v ia^ v iD p iD p iD p iD v h v 
g 2 s hyV u v p F a ^F*h v 



of 

Of 
Of 
Of 
O? 

of 



hy{iD) 2 (iD) 2 hy 
hyiD p ,iD V iD^iD V hy 

hyia^iD^iD^iDfhy 

hyiO^iiDfiD^iDvhy 

h v ia> iV iD li iD p iD v iD p h v 
g 2 s h v F a ^F a uu hy 



(4) 



Below these operators will be denoted collectively by O^ . 

Another 11 operators contain the EOM piece (ivD)h v and are not considered 
any further, since at fixed order l/jrig their matrix elements vanish. 



2 



In addition to the local operators, there are the time-ordered products of the 
lower dimensional operators and . Their contributions at 0(1/ itIq) and 
0(1/toq) are written generically as 



and 



-r(i2) 

ij 

r (lll) 

ijk 



(11) 



(l--Jy)tT 



iT 



n W r0 (2) 



0(1) (D 



z = 1...3. 



i,j = 1...3, i<j 



/n(l) n (i) ro^ 1 ) 



.7 = 1... 7 
i, j,fc = 1. ..3, 



i < j < k . 



(5) 

(6) 
(7) 



The symmetry factor Sijk equals 1, 1/2 or 1/6, if no, two or all inserted operators 
are identical. There is a total of 31 time ordered products but not all of them are 
physical. The reason is that time ordered products, which contain at least one 
operator vanishing by the EOM, can contract to local operators if the (ivD)h v 
term acts on an internal heavy quark line. This is illustrated graphically as 
follows: 

T[[...0 1 {ivD)h v ][h v 2 ...]] oc T[...Oi -0 2 ...] 



Oi 2 Oi ■ 2 

Below such operator relations will be denoted contraction identities (CI). 
They can be derived by manipulating the generating functional of Greensfunc- 
tions in the presence of the relevant operators in much the same way, as one 
usually establishes the validity of the EOM for local operators. As a concrete 
example consider the triple insertion 

T^ 1] =~T [hv(tD) 2 h Vl h v {ivD) 2 h v ~h v {iD) 2 h v ] . (8) 

Naively acting with the EOM operator on the heavy quark propagators to the 
left and the right, one would expect that 

7i ( i3 n) =h v {iDf{iDfh v + ... . (9) 
On the other hand an exact derivation of the contraction identity leads to 

= -iT [h v (ivD)(iD) 2 h v ,h v (iD) 2 h v ] - iT [h v (iD) 2 (ivD)h v ,h v (iD) 2 h v ] 

(10) 

-h v (iD) 2 (iD) 2 h v . 

The apparent contradiction in the sign of the local operator is resolved by the 
fact that the double insertions on the right hand side themselves contain a local 
contribution proportional to h v (iD) 2 (iD) 2 h v : 

iT[h v (ivD)(iD) 2 h v ,h v (iD) 2 h v ] = -h v (iD) 2 (iD) 2 h v + [EOM] (11) 
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Here the symbol [EOM] collects time ordered products that vanish by the EOM 
and can be neclected. The CI allow us to remove all time ordered products 
from the operator basis that contain EOM operator components. That way we 
are left with 4 triple insertions and 4 double insertions collectively denoted as 
q~(m) an d 7~(i2) anc j composed only of the lower order physical operators 0^\ 

(2) 

and 0\ 2- In the next section we will show explicitly that a consistent imple- 
mentation of CI into the renormalization procedure is necessary to maintain 
gauge independence of the renormalized lagrangian. 



3 Anomalous dimensions 

To treat CI correctly during renormalization we start with the calculation of 
the anomalous dimensions corresponding to the full operator basis. This basis 
contains all local operators of dimension 7 allowed by Lorentz invariance and 
the symmetries of HQET as well as all time ordered products which can be con- 
structed from lower order operators disregarding any operator identities (such 
as CI). The anomalous dimension matrix corresponding to the full basis is then 
projected onto the physical basis with all CI applied and EOM operators re- 
moved. The physical basis is minimal in the sense that it contains only linearly 
independent operators. 

To clarify this procedure let us for a moment forget about time ordered prod- 
ucts and consider renormalization of an overcomplete operator basis consisting 
of physical operators Oi and operators Oj which linearely depend on the Of. 

6i = Y,CikO k (12) 

k 

The full basis renormalizes as 

Oi = Z\f Of re + J2 Z^O b 3 are (13) 

j 3 
®i = Yl ^ + Yl ®3 are ' ( 14 ) 
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wherein we use ([L2]) to eleminate the bare operators Oj e in favor of the physical 
ones: 

Oi = £ Zlfof™ + £ Z^c jk O b k are = Z^Of™ (15) 

3 jk 3 

6 % = Z^Of™ + £ zlfc jk O h k a ™ = Z\f0f™ (16) 

3 jk j 

are effective renormalization constants which describe the renormalization 

l 3 

of the physical operator basis Oi among themselves and defines their anomalous 
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~ (2) 

dimensions. The Z>j express the renormalization counterterms of the redun- 
dant operators in terms of the physical ones. They are not linearly independent 
since the renormalized operators on the left hand side of (Jl^Jl^) should also 
fulfill flU), which in turn leads to the consistency condition 

Y,CikZ$ = zg\ (17) 
k 

To one loop order and in the MS scheme, the anomalous dimensions are the 
negative of the pole parts of the renormalization constants: 

Z^ ] =5 l0 -{-)- llJ+ O{{-f) (18) 

Subsequently (|l7|) leads to relations among the anomalous dimensions of the 
physical and unphysical operators, which provides a powerful consistency check 
of the calculation. 

The discussion above applies to our operator basis, if we identify the Oi 
with time ordered products which are redundant after application of the CI and 
the Oi with the residing physical operators. In our case the removal of the 
redundant operators prooves important to ensure gauge independence of the 
physical operator basis. In the following we demonstrate this explicitly in the 
case of the renormalized triple insertion T^H of the kinetic energy operator 
O^. Before applying the CI, the renormalized triple insertion expressed in 
terms of bare operators reads 



2 111 — \~I 



7r e 



C F (2 + T [h v (iD) 2 h v , h v {ivD) 2 h v , h v (iD) 2 h v ] bare 

+ {^C A + C F (™+Z))iT [h v (iD) 2 (ivD)h v X{iD) 2 h v } hare 
+ (i<7 A + C F (y +0)iT [h v (ivD){iD) 2 h v X{iD) 2 h v } hare 
+ {\C A + C F { 1 4 +0) h v (iD) 2 (iD) 2 h h v are 



ig v 3 



(19) 



£ = 1 — £ is the usual gauge parameter and the ellipses denote terms not relevant 
for our discussion. The first term on the right hand side results from local 
renormalization of the kinetic energy operator 

h v (iD) 2 h v = h v {iD) 2 h b v are - (-)-C F (2 + Z)h v (ivD) 2 h b v are . (20) 

The fact that the kinetic energy operator is not renormalized multiplicatively 
is a consequence of the celebrated reparametrization invariance of HQET |l(| . 
However, the kinetic energy operator requires renormalization by the EOM op- 
erator Og 1 ^ , which in turn finds it way into renormalization of the triple insertion 

r (m) 
-'in • 
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The appearance of the operators h v (iD) 2 (ivD)h v and h v (ivD)(iD) 2 h v in the 
double time ordered products on the right hand side of ( |l9| ) derives from the fact 
that the double insertion of the kinetic energy operator needs renormalization 
by a local operator of dimension 6: 

% -T[h v {iDfh v ~h v {iDfh v ) =-fyl(l-C A +C F (™+0 



h v {ivD)(iDyh b v are + h v (Wy(ivD)h% 



bare 



(21) 



If we now replace the bare triple insertion T^ 1 ^ on the right hand side 
of (|l^) with the help of CI (|l0|), the gauge dependencies cancel and we are 
left with a gauge independent result. We stress that the latter property is not 
spoiled by additional application of CI (11) for the double insertions since their 
coefficients are gauge independent already. The result is that apart from the 
above example similar gauge cancellations occur in the case of all time ordered 
products which are related to physical operators by CI. To summarize, as long as 
one works with an overcomplete operator basis in which some of the operators 
are not linearely independent, physical operators may have gauge dependent 
coefficients. However, if one finally projects onto to the physical operator basis, 
these gauge artefacts must drop out. 

We now come to the presentation of our results. The anomalous dimensions 
corresponding to the physical operator basis can be cast into block diagonal 
form 



0(3) 



0(3) 
f (12) 
f-(lll) 



-(3) 

7i 

-(12) 

Ti 

-(in) 



f(12) 



-(12) 
Tnl 
-(HI) 
Tnl 



-f(Hl) 

o \ 





Tnl ' 



(22) 



The first column of this matrix is the result of our calculation and is presented in 
appendix |X| Weinbergs theorem |ll| guarantees, that only local operators are 
needed as renormalization counterterms. However, apart from mixing with local 
operators the triple insertions T^u also require nonlocal renormalization by 



( 121 

the double insertions . This reflects the fact that two operator components 
of the triple insertions arc themselves rcnormalized by a local counterterm of 
dimension 6 (see equation (plf)). 



4 Aspects of the calculation 

To calculate the anomalous dimensions of an operator basis of definite dimen- 
sionality one generally has to insert every operator into potentially UV divergent 
1PI Greenfunctions and extract the pole piece. In the background field method 



G 



Figure 1 : Examples of abelian and nonabelian diagrammatic contributions to 
the renormalization of the triple, double and local insertion. The square, blob 
and triangle represent operators of dimensions 5, 6 and 7, respectively. 



Jl2[ |l3| ], gauge invariance is maintained explicitly, which allows one to get the 
full result from the calculation of a subset of Feyman diagrams with a speci- 
fied number of external background fields. This number is fixed by the highest 
power of gluon field strength tensors appearing in the basis, because such op- 
erators do not contribute at tree level to 1PI Greensfunctions with a smaller 

(3) 

number of external background fields. In our case the local operators C 12/ / 13 
are bilinear in the field strength. This in turn requires the calculation of all one 
loop diagrams with one incoming and outgoing heavy quark and two external 
background fields to get the counterterm contributions of all operators. Figure [I] 
shows some diagrammatic examples. The calculation has been performed with 
the help of algebraic manipulations provided by the package FORM |14j . 

In the first step every integrand corresponding to a specific diagram is ex- 
panded in a Taylor series with respect to its external momenta, i.e. the heavy 
quark momentum and two background field momenta. This procedure is jus- 
tified since we are only interested in UV divergences. The latter are easily 
extracted at this order of the expansion where the power of external momenta 
agrees with the UV degree of divergence of the diagrams. They appear in the 
pole parts of simple one loop tensor integrals of generic structure, 

d d k k^ ■■■ 

(2w) d (vjfc)"»(jfc2)("+4-m)/2 ' ( ' 

which can be related recursively to two basis integrals 7o,o and Iq,2 via Pasarino- 
Veltman reduction. In the calculation integrals with up to 10 indices are needed. 

All diagrammatic contributions to a specific operator insertion are then 
summed up and expressed in terms of tree level contributions of the local oper- 
ators Of ] : 

i=0 

Here A% denotes collectively local operators or time ordered products of the full 
operator basis and (. . . )\pi indicates insertion of the given operator in a 1PI 
Greensfunction to the loop order specified by the superscript. 
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In this step explicit gauge invariance in the background field ensures that 
no gauge variant operators are needed in the series on the right hand side of 
(p4[). This property provides a powerful consistency check of the calculation, 
since the tree level contributions contain up to 30 terms which must match in 
such a way that all gauge variant terms cancel. 

In (|24|) the coefficients 7^ are easily identified as the entries of the anomalous 
dimension matrix 7, , 7; ( ' 12 ^ and 7 ; ( ' 111 ' ) corresponding to mixing with the local 
operators 0\ . In the case of time ordered products the anomalous dimensions 
of their operator components contribute just additively and can therefore be 
derived from the anomalous dimensions of the operator bases appearing at lower 
0(l/m Q ). 

However, up to now the anomalous dimensions correspond to the full opera- 
tor basis and some artificial gauge dependencies of physical operators show up. 
Therefore in the last step of the calculation the projection onto the physical op- 
erator basis has to be performed. This is where unphysical gauge artefacts drop 
out and one is left with the gauge independent anomalous dimension matrix 
of the physical operator basis. These nontrivial cancellations provide another 
powerful cross-check of our results. 



5 Renormalization group logarithms 

With the one loop anomalous dimensions and the tree level matching coefficients 
we are now in the position to solve the renormalization group (RG) equation 
for the Wilson coefficients of the effective lagrangian at 0(1/ m^): 

CW(v)+^ T CW( t i)=0 (25) 



din /i 



Here C (3) (n) denotes the coefficents of the physical operators (<5 (3) , f {12) , f < m )). 
However, an analytical diagonalization of j( 3 > seems to be difficult. Instead 
we restrict ourselves to the calculation of the first logarithmic correction oc 
a s ln(/z/m<Q) in the coefficients C^ 3 \)i). The exact solution of fl25|) reads 



\a s (m Q )J 



(26) 



with the one loop running coupling 



«.(**) =1 _ 2 ^(o) ( ££M )M JL ) (27) 



a s (m Q ) tt tuq' 

where = (33-2n/)/12 in the presence of nf light flavours. Expanding ( p6| ) 
to first order in the strong coupling we get 

C^(„) = C^(m Q ) - (^) H J^W^C^H mQ ) + 0((^n (28) 

TT TUQ TT 
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With our result 7( 3 ) T and the tree level matching coefficients (7' 3 '(mq) the 
Wilson coefficients are easily calculated. The result is shown in the table below. 





tree level 


coefficient of (a s (/i)/7r) ln(/i/mg) 


1 


2 


25/3 C A - 23/3 C F 


2 


-1 


-1/2 Ca 


3 


-1 


4 C A + 8/3 C F 


4 


1 


-17/6C A 


5 


-2 


5/3CU + 8/3CV 


6 


1 


-c A 


7 


1 


-13/6 C A - 8/3 C F 


8 


1 


-Ca 


9 


1 


-4CU 


10 


-1 


9/2 C A 


11 


-1 


9/2 C A 


12 





1/12 C a 


13 





-1/3 Ca 



We only present the coefficients of local operators, since the coefficients of the 
time ordered products are the products of the coefficients of their operator 
components. We note that there is a subtlety in the determination of the tree 
level coefficients C^ {mo). Again the CI forbid their naive extraction from the 
tree level lagrangian at C(l/mg) 

£( 3 ) = h v iD^(ivD) 2 iD tl h v - h v {ivD) A h v - h v ia^ 'iD^ivDfiD^K. (29) 

The reason is that some of the time ordered products with nonvanishing tree 
level coefficients are related to physical operators by CI. Since we are studying 
the RG flow of the physical basis, we first have to reduce the tree level matrix 
elements of the full operator basis with the CI to the ones of the physical basis 
and then apply the RG controlled by the anomalous dimensions of the physical 
basis. 

For example, the triple insertion contributes at tree level with coeffi- 

cient — 1. However, after application of CI (|Tj 11) this coefficient is attributed 
to the local operator = h v (iD) 2 (iD) 2 h v which does not appear in the tree 
level lagrangian. The same occurs for all time ordered products which con- 
tribute at tree level and in addition are affected by a CI. This explains the 
nonvanishing tree level contributions of coefficients of operators not appearing 



(3) 

in the tree level lagrangian. Note that the operators C 12 /i3 
CI contributions because of their color structure. 



are protected from 
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As an phenomenologic application we consider the expansion of the physical 
mass Mjj of a heavy meson in the HQET 



M H = m Q + A - (Ant (OK ' (30) 

with yCi„i given by the powercorrections on the right hand side of (Q) and (• • • } 
denoting the expectation value between heavy meson states in the heavy mass 
limit. Expanding up to 0(1/itiq) we get 

M H = m Q +A- — !— M$ - , n ,, M< 2) - - 1 .. Affi (31) 
y 2m Q ff (2m Q ) 2 ff (2m Q ) 3 ff v 7 



where 



=Y,C1 1 \^)(0^) (32) 



i=l 
2 



< 2) =E c i (9) o*x°* (a) > + E V^Vx^) (33) 

< 3) = Ecf(/x)(0| 3) ) + £ C«(/x)Cf ( M )(7:f )) (34) 

+ E of>(M)^ 1, a.)cl 1 V)<^ u, >. 

i<i<fe=l 

Since the coefficients C^ 2 ^ (fi) are allready known from lower order calculations, 
with our result for the coefficients C\ (//) all short distance contributions up to 
0(1/ (2?tiq) 3 ) are known. However the problem is the nonperturbative input, i.e. 
the hadronic matrix elements on the right hand side of ([52jj3^]34|) which cannot 
be calculated from first principles. The hadronic matrixelements at 0(1/itiq) 
are parametrized by the well known parameters Ai and A2. At C(l/mg) 5 
parameters are needed and at 0(1/toq) a further proliferation of parameters 
takes place: here 21 parameters are needed. 



6 Conclusions 

We have presented a new calculation of the anomalous dimensions of the oper- 
ator basis appearing at 0(1/ nig) of the HQET lagrangian. Local operators as 
well as all time ordered products - composed of lower dimension operators - of 
dimension 7 contribute to the operator basis at this order. We have shown that 
there exist non trivial relations among some time ordered products and local 
operators, so called contraction identities, which have to be properly taken into 
account during renormalization. Otherwise gauge independence of the physical 
operators is violated. 
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The calculation of the one loop anomalous dimensions of the physical op- 
erators was the goal of this paper. Supplemented by the tree level match- 
ing coefficients, it was possible to extract the first logarithmic corrections oc 
a s (/z) In (/z/mQ) to the Wilson coefficents. Of course, to get the non-logarithmic 
contributions at 0(a s ), a one loop matching calculation has to be performed, 
which is beyond the scope of this work. 

In H the short distance coefficients of the 0(l/mg) effective lagrangian have 
been derived by matching the full QCD vertex function with one external gluon 
to HQET in an on-shell renormalization scheme. In this case all loop diagrams 
in the effective theory vanish and the coefficients can be read of directly from 
an expansion of the QCD vertex in powers of the heavy quark momentum and 
the gluon momenta to the appropriate order. It is clear that the coefficnts of 
operators bilinear in the gluon fields cannot be determined by this method. 
Using the equation of motion we have transformed the operator basis of fl into 

(3) (3) (3) 

ours. This way only the coefficients of O y 2 , 6 and O s are determined. The 
other coefficients are related to coefficients which have not been calculated in H . 

(3) 

Unfortunately there is a mismatch in the coefficient of 0% . A possible reason 
for this mismatch may be found in the fact that the tensor structures appearing 
in the expansion of the one loop QCD vertex function at 0(1/toq) cannot be 

(3) 

uniquely attributed to tree level insertions of the local operators 0\ with one 

(3) 

external background field. For example, in our conventions the operators 0\ 

(3) 

and 0\ contribute at tree level identically to the vertex function with only 
one external gluon. Therefore one should expand the one loop QCD vertex 
function with two external background fields in inverse powers of the large mass 
to perform the matching at C(1/wq) correctly. 

In the effect of four fermion operators at 0(l/mg) has been studied. 
The upshot of this work was that the famous Darwin operator is renormalizcd 
by a four fermion operator involving two static and two light quarks, even 
if all four fermion operators are removed at the end of the renormalization 
procedure. Similiar effects cannot be excluded at 0(1/itiq), but in this case the 
higher dimension of the operators make life more complicated. Nevertheless it is 
not only a calculational challenge to include all possible four fermion operators 
into the operatorbasis at C(l/mg) and to investigate their renormalization 
properties 

Another topic which will be covered in a sequel to this publication fl6|j , is 
the reparametrization invariance of HQET. This new symmetry of the HQET 
lagrangian predicts relations among Wilson coefficients of operators appearing 
at different 0(1/itiq). With our results it should be possible to check the validity 
of these relations. This in turn may help to draw a decision between several 
concepts of this symmetry appearing in the literature, all of them unfortunately 
agreeing in their predictions up to 0(l/mQ). 
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A Anomalous dimension matrices 



= Ca a i\ +Cf7; the entries in the first column 



In the decomposition 7 ; 
of (§§) are: 



11 

12 





11 

24 


11 

24 


-4 





1 


-1 


-3 

7 
2 


5 
6 
5 
6 


5 
6 
4 
3 


5 
3 

13 
6 
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7 


1 


1 


1 


1 


6 


8 


12 


8 


3 





7 


1 


5 


1 




8 


2 


24 


3 


1 


5 


-1 


5 


1 


2 


24 




24 


12 





5 


1 


7 


1 




24 


2 


8 


3 


1 


1 

4 





1 

8 





1 


1 


1 


1 


11 


2 


24 


4 


6 


24 


1 


1 


1 


1 


11 


2 


6 


4 


24 


24 


















n 
u 


u 




11 

288 




72 




o 




u 




2 
9 


n 


u 




1 

12 




7 

" 18 


o 


o 




1 

24 




11 

36 


1 

3 


1 

3 












1 
4 


5 

12 












1 

12 


1 

12 












5 

12 


1 
4 












1 
2 


1 
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